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Abstract 

We show that the motion on the n-dimensional ellipsoid is complete integrable 
by exihibiting n integrals in involution. The system is separable at classical and 
quantum level, the separation of classical variables being realized by the inverse of 
the momentum map. This system is a generic one in a new class of n-dimensional 
complete integrable Hamiltonians defined by an arbitrary function f(q,p) invertible 
with respect to momentum p and rational in the coordinate q. 

Complete integrable Hamiltonian systems in the Liouville sense are a main sub- 
ject of interest in the last decades both at classical as well as quantum level. Classical 
examples of such systems are the geodesic flow on the triaxial ellipsoid, Neumannn's 
dynamical system and the integrable cases of the rigid body motion. The number of 
interesting examples increased after the seminal papers by Lax Q and by Olshanet- 
ski and Peremolov especially in connection with the inverse scattering method 
Q\ , and the relation with classical simple Lie algebras. 

The aim of this paper is the study of the geodesic motion on the n-dimensional 
ellipsoid, which is a direct generalization of the 2-dimensional case studied by Jacobi 
|| , since it seems that no solution is known for this system in the case n > 2 . We 
obtain a complete description of the problem at the classical level by finding n 
prime integrals in involution, the separation of variables, the explicit solution of the 
Hamilton- Jacobi equation and the equations of geodesies. We also show that the 
Schrodinger equation separates. This system is to a great extent universal among 
the various integrable systems and generates a new class of completely integrable 
models. 

The Lagrangean for a particle of unit mass constrained to move on the n- 
dimensional ellipsoid 

fl + f2+.. . + ^±1 = 1 (1) 
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The preceding equations define a constrained system and the model can be for- 
mulated in terms of constrained dynamical variables with Dirac brackets, or in an 
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unconstrained form with canonical Poisson brackets. In the following we will use a 
third way, the classical one, which consists in reducing the Lagrangean by eliminat- 
ing one degree of freedom || using the equation of constraint (1). The simplest way 
to do that would be to resolve Eq.(l) with respect to the last coordinate and use it 
in the free Lagrangean, Eq.(2), but the drawback is that we obtain a non-diagonal 
metric. We follow here the Jacobi idea which was to find a clever parametrization 
such that the coresponding metric should be diagonal ||, |6|. 
For what follows it is useful to define two polynomials 

n+1 n 

P(x) = I]> - tH) Q(x) = Y[(x-ih) (3) 

i=i i=i 

where a, and Ui, i = 1, . . . ,n are the positive numbers entering the parametriza- 
tion of the ellipsoid and the ellipsoidal coordinates, respectively. The orthogonal 
parametrization of the quadric (1) is given by @ 



where P'(aj) = dP{x)/d 

x\x=dj •> and the ellipsoidal coordinates txi, . . . ,u n satisfy 
a± < u± < a>2 < ■ ■ ■ < u n < a n +i- Using this parametrisation in Eq.(2) we find 
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where the (diagonal) metric is given by ga = Ui Q' (u^/ P(ui) ,i = 1,2, ... ,n, and 
Q'{ui) = dQ(x)/dx\ x=Ui . Defining as usual the generalized momenta by pi = 
d£/diii and using the Legendre transform we find the Hamiltonian of the problem 



i=n i=n 
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where g n = P{u,i)/ui Q'{u; L ). We define now the symmetric functions of the polyno- 
mials QV\x) = Q(x)/(x — Uj) 



n-l 

Q«\x) = J2x k si;l k _ 1 , j = l,2...,n (7) 



fe=0 

where Sq^ = 1, = — {u\ + . . . + Uj-i + Uj + \ + . . . + u n ), etc. The upper index 

means that the coordinate Uj does not enter the symmetric sum S^',k = I,..., n-l. 
We define the following functions 

Hk = Y = E S ti Pi k = l,2...,n (8) 

i=i i=i Ui ^ {Ui) 

where Hi differs from Sj by a numerical factor. A careful inspection of the Eqs.(8) 
shows that for each degree of freedom the contribution to the Hamiltonian 
is given by the product of two different factors. The first one depends on the 
"Vandermonde" structure f± = S^^/Q' \ui) and the second one fa = (P(u{)/ui)pf 
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depends on the "singularities", i.e. the hyperellitic curve defined by the parameters 

Let g(p,u) = 7i(p, u) be an arbitrary function depending on the canonical vari- 
ables p and u which is invertible with respect to the momentum p. As we will see 
later the invertibility condition is necessary for the separation of variables in the 
Hamilton- Jacobi equation. In particular we may suppose that TC(p,u) is an one- 
dimensional Hamiltonian. For each n 6 N we define an n-dimensional integrable 
model by giving n integrals in involution 

n 

^-k(p,u) = J27W^\9{Pi,Ui), k = l,2,...,n (8') 

i=1 V \ u i) 

Our main result is contained in the following proposition. 

Proposition. Lei M 2n ~ T*(R n ) be the canonically symplectic phase space of 
the dynamical system defined by the Hamilton function Eq.(6). Then 

i) the functions H i: i = 1, 2, . . . , n are in involution 

{H i ,H j } = 0, i,j = 1,2,..., n 

ii) the momentum map is given by 

£:M 2n ^R n :M h = {(u l , Pi ) : H t = K, i = 1, 2, . . . , n}, /i,Gl 

then £ _1 (Mh) realizes the separation of variables giving an explicit factorisation of 

Liouville's tori into one- dimensional ovals 

Hi) canonical equations are integrable by quadratures 

In the following we sketch a proof of the above proposition. As it will be easily 
seen the same proof is also true for the Hamiltonians defined in Eqs.(8'). 
Proof. 



i) We calculate the Poisson bracket 



{H k ,H l } = Y J 
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The last step was possible because the symmetric functions Si and depend 
on all u\, U2, ■ ■ ■ , u n , but Uj. Looking at the last expression it is easily seen that the 
partial derivative with respect to Uj vanishes for i = j. For i / j we have to show 
that 

r, q(j) q(i) q(j) 
° 1-1° k-1 ~ °k-l°l-l 



duj 
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but this is a direct consequence of the following identities 

where the upper index means that the corresponding expression does not 

depend on both U{ and uj. In this way we have shown that {H^, Hi} = 0. 

ii) The surface Mh = {(ui,pi) : Hi = hi, i = 1,2, ... ,n}, hi € R is given by 
the system of equations 
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Pi=h k , k = l,2...,n 
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and resolving it with respect to pi is equivalent to the calculation of the following 
determinant 
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which is equal to the Vandermonde determinant, i.e. 

V n = Yi ( U 3 ~ Ui ^> 
l<i<j<n 

(i) 

Let V^_i be the determinant obtained by removing the jth column and the last 
row in V n and W n j the determinant obtained by replacing the jth column of V n by 
(hi, . . . , hnf . It is easily seen that 

v^i= n («•-«*) 

l<fc<!<n 



i.e. Vn_i is the Vandermonde determinant of the variables ui, . . . ,u n , but uj. We 
have the identities 
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Using these identities it is easily seen that £ 1 (M\ l ) is equivalent to the relations 

n-l 



P(ui)p: 



u. 



= ^2 K-kUi i = l,2... 



n 



(10) 



fc=0 



which shows that £ 1 (M\ i ) is a n-dimensional submanifold of M 2n and more impor- 
tant pf are functions which depend only on the variable Ui. The last relation shows 
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that the application 6~ l (M\ i ) realizes the separation of variables for the geodesic 
motion on the ellipsoid. 

For the Hamiltonians given by Eqs.(8') the preceding equations have the form 



n-l 



g{pi,ui) = K-kUi i = l,2...,n (10') 



k=0 

The relations (10-10') have the classical form || 

(p(xi,pi, hi, . . . , h n ) = i = \,2...,n 

which is an explicit factorization of Liouville's tori into one-dimensional ovals. This 
is important for the quantization problem. 

With the notation R(u) = ^2^=0 h-n-k uk the above relations can be written 



UiR(ui) . 
Pi = e i\ ni \ > i = !,••■, n 

Pi = g~ 1 (R(u i )), i = l,...,n 

where ej = ±1 and g _1 is the inverse of the relation (10') with respect to the 
momentum p. The last relations allow immediately to resolve the Hamilton-Jacobi 
equation because in this case the action splits into a sum of terms 

S(h,ui, ...,«„) = Si(h,ui) + . . . + S n (h,u n ) 

each of them satisfying an ordinary differential equation. Only the solution of the 
Hamilton-Jacobi equation for the geodesic motion on the ellipsoid will be presented 
the other more general case being similar. 



S(h,m, . . . ,u n ) = V €i / J dw 

Hi) The above formulae allows us to choose new canonical variables as follows 
Qi = Hi, Qk = Hk, k = 2, . . . , n and the corresponding variables Vi, i = 1, . . . ,n. 
The Hamilton equations take the form 

Qi = 0, i = 1, ... ,n 

Vi = -\, Vi = 0, i = 2,...,n 

and therefore Qi = hi, i = 1, . . . , n and V\ = — t + g%, Vk = gki k = 2, . . . , n, with 
gi,hi £ R, i = 1, . . . ,n. 
Because 

as _ ds _ ei r Ui t n - i+l > 2 

% ~ ~~dQ~~ ~dhi ~ 2 J u o y/P(t)R(t) 1 
we obtain the system 

—tdu + bn = - > e, / — = dt, j = 1, .... n 
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which gives the implicit equations of the geodesies. In this way the integration of 
the Hamilton equations was reduced to quadratures. On the last expressions one 
can see that all the subtleties of the geodesic motion on the n-dimensional ellipsoid 
are encoded by the hyperelliptic curve y 2 = P(x) R(x) whose genus is g = n. 

For quantization we use both the forms (9) and (10) and show first that the 
quantization of Hi is equivalent to the quantization of Eq.(10). It is well known 
that because of the ambiguities concerning the ordering of p and u we must use 
the Laplace-Beltrami operator ||. Its general form is A n = -^Pi(^ g^Pj), i,j = 
1, . . . , n, where g = det(gij) and g^ is the metric tensor. Taking into account that 

■ V v- 

g ii = u i Q'{u i )/P{u i ) = {-l) n - " ' 



V® P( Ui 

n— 1 



after some simplifications the Schrodinger equation generated by the Hamiltonian 
Hi is written in the following form 



(i) 

Since the term does not depend on itj it can be pulled out of the bracket and 
the precedent equation takes the form 



i=l 



Now we make use of the Jacobi identity for the Vandermonde determinant. With 
the above notations the identity is 



rti) „,fc 

3 

and using it in the preceding relation we obtain 

n-1 



8=1 



P(Uj) d P(uj) d*\ u ^ 

m dui \ V m dm l ^ n ~ k 1 



k=0 

which is equivalent with n indepedent equations of the form 

'n-1 \ 



ip( Ui ) d / p( Ui )d^i 



m dui \ \ Ui dm 



+ [J^Cn-kUi )*i = 0, i = l...,n (11) 



vfe=0 



Here c\ = hi and the other are arbitrary. The direct approach, starting from 
Eq.(10), is simpler the problem being one-dimensional and one arrives at the same 
equation, Eq.(ll). It has the advantage that the arbitrary coefficients are iden- 
tified to Cfc = hk, i.e. Cfc are the eigenvalues of the Hamiltonians H^. 
The Eq.(ll) has the general form of a Sturm-Liouville problem 

4(^W)/(*HM*),(*) 
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which has to be resolved on an interval [a, b\. It is well known that its eigenfunctions 
will live in a Hilbert space iff p(x)r{x) > on [a, b]. If p(x) has a continuous 
first derivative and p(x) r(x) a continuous second derivative, then by the following 
coordinate and function transforms 

<P = f (^||) V2 dx, $ = (KuMu)) 1 /* f(u) (12) 

the preceding equation takes the standard form 

d 2 <S> 



where 



+ q{v)<& = A<I> 



fi ((p) v(u) 1/4 

W = — pr, Mv) = P " r(u)) 1 

H{(p) r{u) 

and u = u{ip) is the solution of the Jacobi inverse problem (12). 
In our case, Eq.(ll), the transformation is 



/"( 



u R(u) 



1/2 



and the Schrodinger equation has the form 

where /x(y?) = (^(^(v 7 ))) 1 ^ 4 an d in -R(u) we made the rescaling hk — > h^/hi, k = 
1, . . . ,n. Thus we have obtained that the solving of the Schrodinger equation (11) 
is equivalent to the solving of the motion of one-dimensional particle in a potential 
generated by R(u(ip)). 

For n = 1 Eq.(13) is the equation of the one-dimensional rotator 

with the solution Hf((p) = e tl(f , I £ Z, etc. In all the other cases we have to make 
use of the theory of hyperelliptic curves, ^-functions and/or hyperelliptic Abelian 
functions in order to obtain explicit solutions. This problem will be treated else- 
where. What is remarkable is that the solving of the classical problem, or the solving 
of the associated Schrodinger equation leads to the use of the same mathematical 
formalism, ^-functions or hyperelliptic Kleinian functions. 

However there is a simpler alternative to Eq.(13), namely by the change of 
variable 

* = f fLdt 



p(t) 



the Eq.(ll) takes the form 

d 2 ^ 



+ l^2h n ^ k u(ip) k \ * 



\k=0 
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which is the equation of a particle moving in a potential generated by the integrals 
in involution. 

As concerns the quantization of the Hamiltonian H i it depends on its explicit 
form and we do not pursue it here. 

From the proof of our results it follows that the hyperelliptic curve was only a 
tool in obtaining the separation of variables Eq.(lO); in fact the separation was a 
direct consequence of the properties of the Vandermonde determinant. 

In the following we exibit a few examples of new n-dimensional integrable mod- 
els. Two models which show that the dimension n of the system has no direct 
connection with the number of zeros and/or poles of the function g(p,u) could 
be: g(p,u) = (sinu/u) p 2 and g(p,u) = tgue ap , the first example being a func- 
tion with a denumerable number of zeros and the second one with a denumer- 
able number of poles and zeros, in both cases the hyperelliptic curve being of infi- 
nite genus. Other examples are deduced for example from the many-body elliptic 
Calogero-Moser [flO| ], or the elliptic Ruijenaars models pj| . Starting with the corre- 
sponding one-dimensional Hamiltonians Hcm(p,u) = f(p,u) = p 2 /2 + u 2 p T (u) and 
fn(p,u) = cosh(ap)^Jl — 2 {au) 2 p T {u) respectively, where p T (u) is the Weierstrass 
function, we obtain n-dimensional models. 

In conclusion we discovered in this paper a new class of complete integrable sys- 
tems which allows to uncover the origin of their integrability or solvability property. 
We have shown that there is a simple and general mechanism allowing us to con- 
struct complete integrable Hamiltonian systems with an arbitrary number of degree 
of freedom, and for all these systems the separation of classical variables is given by 
the inverse of the momentum map. 
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